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I.  Status  of  Effort 

In  this  project  the  dissimilarity  between  streamwise  velocity  and  temperature  in  developed  or  non- 
perturbed  (where  ’developed  condition’  means  a  region  downstream  of  a  point  beyond  which  the  flow 
field’s  behavior  is  streamwise  invariant)  and  perturbed  turbulent  channel  and  plane  Couette  flows,  was 
addressed  using  direct  numerical  simulations,  DNS.  The  main  goal  of  this  study  addressing  the  dis¬ 
similarity  between  streamwise  velocity  and  temperature,  was  to  obtain  insights  that  can  aid  turbulent 
heat-transfer  modeling  based  on  Reynolds  stress  (Note  that  throughout  this  report  the  expression  mean- 
field  dissimilarity  will  be  used  for  the  difference  ( U  —  0)  in  dimensionless  form,  where  U  and  0  are  the 
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mean  streamwise  velocity  and  temperature,  respectively;  dissimilarity  of  streamwise  velocity  and  tem¬ 
perature  fluctuations  for  the  dimensionless  difference  ( u'  —  9')',  streamwise  turbulent  flux  dissimilarity 
for  the  dimensionless  difference  ( (v'v')  —  (• Q'u '));  and  wall-normal  turbulent  flux  dissimilarity  for  the 
dimensionless  difference  {{u'v')  —  ( 6'v '))). 

First  the  natural  dissimilarity  of  streamwise  velocity  and  temperature  fluctuations,  in  non-perturbed 
turbulent  channel  and  plane  Couette  flows,  was  characterized  (Kestin  and  Richardson,  1963;  Fulachier 
and  Dumas,  1976;  Antonia  et.  al,  1987).  For  both  of  the  flow  configurations,  a  Reynolds  number  based 
on  the  friction  velocity  of  about  150  was  used  (Pasinato,  H.D.,  Velocity  and  temperature  dissimilarity  in 
fully  developed  turbulent  channel  and  plane  Couette  Flows,  Int.  J.  Heat  F.  Flow,  32,  pp.  1 1-25,  201 1). 

In  the  second  paid  the  dissimilarity  between  the  mean  turbulent  fluxes  of  streamwise  velocity  and 
temperature  in  perturbed  turbulent  flows  was  characterized  (Kong,  Choi,  and  Lee,  2001).  Numerical 
experiments  for  turbulent  channel  and  plane  Couette  flows  were  implemented,  with  analogous  mean 
streamwise  velocities  and  temperatures  in  developed  conditions,  which  were  then  perturbed  in  different 
ways  in  a  second  DNS.  To  provide  inflow  boundary  conditions  for  the  perturbed-developing  flow,  a 
parallel  DNS  with  identical  numerical  resolution  was  performed.  The  perturbations  used  in  the  study 
were  (a)  local  blowing  or  suction  from  a  narrow  span-wise  slot  at  the  walls,  (b)  local  streamwise  pressure- 
gradient  steps,  adverse  and  favorable,  in  a  naixow  span-wise  volume  in  the  buffer  region,  (c)  local  wall- 
normal  pressure-gradient  steps  in  a  naixow  span-wise  volume  in  the  buffer  region,  and  (d)  local  wall 
temperature  steps.  The  main  idea  was  to  study  the  disruption  of  mean-field  similarity  due  to  turbulence 
in  perturbed  turbulent  flows  and  to  study  the  main  causes  of  dissimilarity  of  the  turbulent  fluxes  for 
streamwise  momentum  and  heat  (Pasinato,  H.D.,  Streamwise  Velocity  and  Temperature  Dissimilarity  in 
Perturbed  Channel  and  Plane  Couette  flows,  to  be  submitted). 

II.  Accomplishments  and  New  Findings 

=>  It  was  found  that  the  natural  dissimilarities  were  more  associated  with  those  events  in  (Q 2  +  Qf), 
according  to  the  quadrant  analysis  technique  (at  about  60%  )  for  both  flow  configurations.  In 
other  words,  the  same  types  of  events  responsible  for  wall-normal  turbulent  fluxes  (streamwise 
momentum  and  heat)  are  also  responsible  in  paid  for  the  dissimilarities. 

=>  In  the  frequency  domain  it  was  shown  that  the  dissimilarity  is  caused  by  a  shift  toward  higher 
frequencies  as  the  distance  from  the  wall  increases  for  the  9'  spectrum  in  comparison  with  the  v! 
spectrum. 


ii 


=>  At  the  wall-layer  the  O'  spectrum  presented  higher  frequencies  than  u!  and  lower  frequencies  than 
p'  (and  also  than  dp' /dx  spectrum  ).  These  differences  generate  a  more  significant  correlation  of 
O'  with  dp' /dx  than  u!  with  dp' /dx.  This  correlation  difference  is  the  main  cause  of  the  dissimi¬ 
larities  in  non-perturbed  turbulent  flows;  or  in  the  decorrelation  of  v!  and  O'  from  the  wall  toward 
the  central  region. 

=>  For  perturbed  turbulent  flow  the  dissimilarity  for  both  flow  configurations  were  analogous. 

=>  The  contribution  to  the  mean- field  dissimilarity  by  turbulent  fluxes  dissimilarity  was  a  minor  part 
in  comparison  with  the  pressure  gradient  contribution. 

=>  The  contribution  to  mean-field  dissimilarity  by  the  streamwise  turbulent  fluxes,  ( (u'u')  —  ( u'O ')), 
mainly  occurs  in  the  perturbation  region,  while  those  due  to  wall-normal  fluxes,  ((v'u')  —  (v'O')), 
occurs  mostly  downstream  of  the  perturbation  region. 

=>  The  wall-normal  turbulent  fluxes  make  a  greater  contribution  to  mean-field  dissimilarity  than  the 
streamwise  fluxes. 

= >  The  leading  dissimilarity  contributor  to  the  wall-normal  turbulent  fluxes  dissimilarity  was  (v'v')d<& / dy, 
where  T  is  the  difference  (U  —  0). 

=>  A  simple  transformation  of  the  Reynolds  stresses  in  the  streamwise  and  wall-normal  direction, 

(■ u'u '),  (v'v'),  gives  a  reasonable  a  prior  prediction  of  the  turbulent  heat  fluxes  ( O'u '),  (O'v1). 

III.  Natural  dissimilarity  in  non-perturbed  turbulent  flows 

(III.  a)  Definition  of  a  new  variable  to  measure  dissimilarity 

In  the  numerical  experiments  the  streamwise  velocity  and  temperature  are  analogous  in  the  mean  values 
(the  Reynolds  averaged  conservation  laws  for  streamwise  velocity  and  temperature,  with  some  restric¬ 
tions  as  constant  fluid  properties,  Pr  =  1,  and  identical  boundary  conditions,  arc  analogous).  Thus 
(j)  =  (U  —  0)  +  (v!  —  O'),  resulting  that  =  <b  +  ft  as  all  variables,  where  the  variance  of  ft  is 

VAR((f>')  =  (</>+)2  =  (u'u')  -  (u'O')  +  (O' O')  -  (u'O').  (1) 

The  local  variance  V  AR(f)  was  used  as  a  measure  of  the  mean  local  dissimilarity  and  the  averaged 
variance  V AR((j)')ave  as  an  averaged  value  in  the  whole  domain. 
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Figure  1:  Percentage  of  V AR(<f>')  for  events  with  different  conditions.  Solid  line,  P(<f>'  2  >  kf+2;  u'v'  <  0); 
-  -  -  -  P(ft  2  >  fc</>+2;dV  <  0);  ...  .  ....  Pff2  >  k(j>+2]v/v'  <  0;P  <  0);  □.□.□.□  , 

P{4>'2  >  k(j>+2\6'v'  <  0;  v'  <  0) ;  >•>•>•>,  P{4>'2  >  k(j)+2\ dp'/dx  >  0).  (a)  Channel  flow;  (b) 
Couette  flow. 

(lll.b)  Characterization  of  a  wail -layer  event  associated  with  dissimilarity  in  non-perturbed  turbulent 
flows 

The  relationship  between  V AR{f)  and  those  characteristic  events  in  the  wall  layer  responsible  for  the 
wall-normal  turbulent  transport  was  studied  using  conditional  probability.  The  association  of  these  events 
with  the  dissimilarity  was  evaluated  by  first  detecting  events  that  satisfy  the  condition  ft2  >  kf+2,  where 
(j)+  is  a  long-term  statistic  from  a  previous  large  sample,  k  is  a  parameter,  and  (j/2  is  the  variance  of  the 
dissimilarity  during  this  event.  Then,  once  an  event  with  this  condition  was  detected,  which  for  simplicity 
was  called  an  important  dissimilarity  event  (IDE),  a  second  (or  even  a  third)  condition  (e.g.,  the  sign  of 
the  wall-normal  velocity  v'  <  0,  etc.)  was  evaluated,  and  its  dissimilarity  V AR(f/) [ was  evaluated. 
A  detection  algorithm,  analogous  to  those  used  in  the  literature  to  detect  burst  or  ejection  events,  was 
used  (Blackwelder  and  Haritonidis,  1983;  Johansson,  Her  and  Haritonidis,  1987;  Luchik  and  Tiederman, 
1987).  In  this  study,  however,  the  goal  was  not  to  detect  these  kinds  of  events  or  evaluate  the  V AR{f) 
associated  with  them.  Rather,  the  goals  were  to  detect  a  high  amplitude  in  <f>'  greater  than  some  value, 
or  an  IDE,  checking  whether  it  was,  for  example,  an  event  in  the  second  or  fourth  quadrant  (Q 2  +  Q 4) 
according  to  the  quadrant  analysis  technique,  for  Reynolds  stress  or  wall-normal  turbulent  heat,  etc.,  and 
to  then  evaluate  its  contribution  to  V AR{f)ioE- 
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VAR  =  (j)' 2  >  kcj)+ 2 


(2) 


where  the  mean,  < f>,  and  rms,  f+,  values  used  in  previous  condition  (2)  were  long-term  statistics,  evalu¬ 
ated  from  a  previous  long  sample,  and  the  wide-hat  symbol  means  a  short  term  mean  value  in  the  T  time 
filtering  interval. 


- — ~ — -  1  rt+T/2 

4>(t,T)2  =  —  /  ftrftdr  (3) 

1  Jt-T/2 

Then,  once  an  event  with  condition  (2)  was  detected,  conditional  probability  with  different  conditions 
was  applied  to  determine  whether  the  IDE  satisfied  a  second  or  a  second  and  third  condition.  Some  of  the 
conditions  included  the  following:  P(ft 2  >  kcj)+2,  u'v'  <  0),  to  determine  whether  these  events  also  be¬ 
longed  to  events  in  the  second  or  fourth  quadrant  (Q2+Q4)  for  Reynolds  stress;  P(ft 2  >  k(f>+2,  u'v'  <  0, 
to  determine  whether  these  events  also  belonged  to  events  in  the  fourth  quadrant  (Q 4)  for  Reynolds 
stress;  or  the  condition  P(ft 2  >  left'  2,  u'v1  <  0,  9'v'  <  0),  to  determine  whether  the  dissimilarity  asso¬ 
ciated  with  events  in  (Q 2  +  Q 4)  for  Reynolds  stress  was  close  to  the  dissimilarity  associated  with  events 
in  (Q 2  +  QA)  for  Reynolds  stress  and  wall-normal  turbulent  heat  at  the  same  time.  Figures  1(a)  and  1(b) 
show  some  results  from  this  algorithms. 

(III.c)  Natural  dissimilarity  in  the  frequency  domain 

Using  the  new  variable  the  distribution  of  energy  of  the  different  spectra  was  evaluated.  Fig.  2(a)  and 
2(b)  show  v! ,  O' ,  //,  and  ft  spectra.  These  figures  show  that  u'  and  O'  spectra  present,  for  both  flow 
configurations,  a  shift  toward  higher  frequencies  as  the  distance  from  the  wall  increase,  being  always 
more  important  for  O'  than  for  u'  (this  behavior  is  clearer  from  the  spectrum  of  their  difference,  ft). 

(Ill.d)  Velocity  and  temperature  pressure-gradient  interaction 

There  is  a  more  significant  correlation  of  O'  with  dp' /dx  than  v!  with  dp' /dx  (Figure  3).  This  correlation 
difference  is  the  main  cause  of  the  dissimilarities  in  non-perturbed  turbulent  flows;  or  in  the  decorrelation 
of  u!  and  O'  from  the  wall  toward  the  central  region. 

IV.  Dissimilarity  in  perturbed  turbulent  flows 

(IV.a)  Scheme  of  numerical  simulations 

To  disrupt  the  mean-field  analogy  that  exists  in  the  developed  conditions  of  the  first  DNS,  different  kinds 
of  perturbations  were  used  in  a  narrow  region  of  the  computational  domain,  or  at  the  boundaries,  of  the 
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Figure  2:  Spectral  density  functions  of  u' ,  9',  <//,  and  //  at  four  positions  from  the  wall  for  (a)  channel  and  (b) 
Couette  flow.  Viscous  region  (top-left);  buffer  region  (top-right);  beginning  of  logarithmic  region  (bottom-left); 
center  region  (bottom-right).  Solid  line,  a  =  u'/u+\  o  •  o  •  o  •  o ,  a  =  9' /9+;  a  =  </>'/</>+; 

a  =  p1  /p+ .  Note  the  different  scales. 
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Figure  3:  —  (</>' dp' /dx).  —  .  —  .  — ,  total;  solid  line,  at  Q4ruiv>y, - ,  at  Q2(u>v>y,  o  •  o  •  o  •  o,  Q4^g 

Q2(givi\.  Non-dimensionalized  with  u\jv.  Note  the  shift  in  the  ordinate,  0.3,  for  the  channel  plot. 


Developed  flow 


Developing  flow 


Figure  4:  Scheme  of  parallel  DNS  for  a  numerical  experiment  of  perturbed  turbulent  flow. 


Figure  5 :  Domain  and  location  of  spanwise  slot  used  for  perturbation.  Solid  vertical  lines  denote  the  slot  location; 
vertical  dotted  lines  denote  sections  where  solution  is  presented,  and  horizontal  dotted  lines,  positions  from  the 
wall  where  solution  is  presented.  (a)Channel;  (b)Couette  flow. 

second  DNS.  These  perturbations  were  (a)  local  blowing  and  suction  from  a  narrow  span-wise  slot  at  the 
wall,  (b)  local  adverse  and  favorable  streamwise  pressure-gradient  steps  at  a  narrow  span-wise  volume 
of  the  buffer  region,  (c)  local  adverse  and  favorable  wall-normal  pressure-gradient  steps  at  a  narrow 
span-wise  volume  of  the  buffer  region,  and  (d)  local  wall  temperature  steps. 

The  following  is  a  list  of  the  different  numerical  experiments  reported  in  publications,  which  present 
the  most  important  dissimilarity  between  turbulent  fluxes. 

1.  CFIB6220:  Perturbed  channel  flow  with  blowing/suction  from  the  lower  wall,  W+  =  220,  v+  = 
±0.60.  Where  W+  is  the  width  of  the  slot  and  v+  the  injected  wall-normal  velocity. 

2.  CFB6590:  Perturbed  plane  Couette  flow  with  blowing/suction  from  the  lower  wall,  W+  =  590, 
v+  =  ±0.60. 

3.  CHABR220:  Perturbed  channel  flow  with  an  adverse/favorable  streamwise  pressure-gradient  step 
at  the  buffer  region  of  the  lower  half,  W+  =  220,  (dP/dx)+  =  0.30  for  adverse  and  —0.30  for 
favorable. 

4.  CFABR590:  Perturbed  plane  Couette  flow  with  an  adverse  pressure  gradient  step  at  the  buffer 
region  of  the  lower  half,  W+  =  590,  {dP/dx)+  =  0.18  for  adverse  and  —0.18  for  favorable. 

(IV.b)  Leading  contributor  terms  to  turbulent  fluxes  dissimilarity 
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Figure  6:  Contribution  to  differences  (< f>'u ')  (lhs)  and  ((f>'v')  (rhs)  by  turbulence  production  for  cases  (a) 
CHB6220,  (b)  CHABR220,  (c)  CFB6590,  and  (d)  CFABR590  at  y+  =  38.  +•  +  •  +  ■  +,  total;  solid  line,  function 
of  mean-field  dissimilarity; - ,  function  of  turbulence  dissimilarity. 
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The  differences  between  (( u'u ')  —  ( Q'u '))  depend  primarily  on  mean  flow,  and  secondarily  on  turbulence, 
with  (u'v')d<& /dy  as  the  leading  term.  Conversely,  the  differences  between  (( u'v ')  —  {Q'v'))  depend 
mostly  on  mean  flow  dissimilarity,  with  (v'v')d<&  /  dy  as  the  leading  term.  In  other  words  mean  fields 
wall-normal  gradients  play  the  most  important  rule  in  the  dissimilarity  between  streamwise  velocity 
and  temperature  mean  turbulent  fluxes.  Figures  6(a),  6(b),  6(c),  6(d)  show  these  results  for  blowing 
and  adverse  pressure  gradient  step  for  both  flow  configurations.  In  those  regions  of  the  flow  where 
d(\> /dy  =  dU/dy  —  dQ/dy  =  0  dissimilarity  is  minimal.  Both  flow  configurations,  channel  and 
Couette  flow,  present  similar  results. 

(IV.c)A  relation  that  holds  for  non-perturbed  and  perturbed  turbulent  flow 


Figure  7:  9rrns/urms  and  (dQ  /  dy)  /  (dU  /  dy)  for  channel  flow  perturbed  with  (a)  blowing  and  (b)  temperature 

step  at  the  wall,  o  •  o  •  o  •  o,  non-perturbed  values.  Solid  line,  on  the  slot; - ,  W+ ;  —  .  —  .  —  5W+ 

downstream. 


For  Pr  =  1  the  following  relation  holds  for  the  wall-layer  for  the  simple  non-perturbed  and  perturbed 
turbulent  flows  used  in  this  study  (Figures  7(a)  and  7(b)): 


^ rms  Urms 

dU  /  dy  <90  /  dy 


(4) 


(IV.  d)  Turbulent  heat  transfer  modeling 

The  following  transformations  of  the  Reynolds  stresses  using  the  wall-normal  gradients  of  streamwise 
velocity  and  temperature  show  a  reasonable  a  priori  predictions  of  the  turbulent  heat  fluxes,  as  it  is  shown 
in  Figures  8(a)  and  8(b),  for  both  flow  configurations. 
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Figure  8:  A  prior  comparison  of  modeled  turbulent  heat  flux  (O' v!)  and  {0'v')+ ,  for  channel  flow  perturbed  with 
(a)  blowing  and  (b)  an  adverse  axial  pressure  gradient  step.  Solid  line,  and  +  ■  +  ■  +  ■  +(second  symbol  modeled), 

on  the  slot; - ,  and  x  •  x  •  x  •  x,  W+ ;  —  .  —  .  —  and  .  5 W+  downstream;  o  •  o  •  o  •  o, 

non-perturbed  values. 
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where  all  variables  are  in  dimensionless  form. 
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